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This paper considers the computational capabilities of different 
mathematical models of magnetic bubble interactions. A specific model 
was studied earlier by R. L. Graham, who showed that there exist 
combinational functions of 11 or more variables that cannot be com- 
puted by this model. This paper extends his results by introducing 
different types of interactions which seem to be practical and enable 
the computation of all combinational functions. The problem of 
efficient computation from the points of view of time and space re- 
quirements and the geometrical requirements imposed by the fact that 
interactions can occur only between physically adjacent locations are 
also examined. Finally, a model in which computations are carried 
out by applying uniform magnetic fields to the entire platelet, with 
individual access limited to locations along the periphery, is presented. 

I. INTRODUCTION 

Cylindrical magnetic domains in certain orthoferrite materials have 
been investigated extensively in recent years. 1-3 These domains, com- 
monly referred to as bubbles, have the property that they can be 
moved within the material by the application of suitable external 
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magnetic fields. The motion of a bubble is also dependent on other 
bubbles in its vicinity. The locations of bubbles can be restricted to 
a finite set of possible positions in an orthoferrite platelet. The pres- 
ence or absence of a bubble at a particular location may be treated as 
representing the values of a binary variable. Thus, magnetic bubbles 
seem to have natural applications in performing memory and logic 
functions. 

In a recent paper, R. L. Graham 4 considered the computational 
capabilities of a particular mathematical model of magnetic bubble 
interaction. The only type of interaction allowed in this model is 
represented by an instruction of the type {x, y) where x and y are 
distinct adjacent locations in the orthoferrite platelet. Application of 
this instruction results in moving the bubble in location x to location y 
if the latter does not contain a bubble prior to the application of the 
instruction. If x does not contain any bubble or y already has a bubble, 
no transfer of bubbles takes place. Two locations, only one of which 
contains a bubble, are used to represent each binary variable and its 
complement. Graham has shown that only a small fraction of all com- 
binational functions of 11 or more variables can be computed by this 
model. 

In this paper, we review Graham's results and extend his model 
by introducing different types of interactions which seem to be prac- 
tical and enable the computation of all combinational functions. We 
then examine the problem of efficient computation from the points of 
view of space and time requirements. The geometrical requirements, 
necessitated by the fact that bubble interactions can occur only be- 
tween physically neighboring locations, are also examined. Finally, 
we consider a model in which computations are carried out by apply- 
ing magnetic fields to the entire platelet and individual access is lim- 
ited to locations along the periphery. 

II: MODELS OF BUBBLE INTERACTIONS 

Let us first consider the model studied by Graham. 4 All locations are 
assumed to be adjacent to one another so that interaction between any 
pair of locations is possible. If S is the set of all possible bubble loca- 
tions and n is the number of such locations, let Q be the set of all 2" 
subsets of S. If a subset X tQ containing bubbles represents the values 
of the variables of a function /, then / is a mapping from Q onto Q. 
The function / is realized by applying a program P to S with bubbles 
in locations contained in X. If the set of resultant bubble locations is 
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X P , then X p = f(X) and X p t Q. The program P consists of a sequence 
of instructions of the form e = (a, b) such that X e = (X — {a})VJ {&}, 
if a e X, b 4 X, and X" = X otherwise. The nondecreasing overlap 
(NDO) theorem due to Shockley and proven in Ref. 4 points out some 
of the limitations of this model of interaction. The NDO theorem is 
repeated below in the interest of completeness. 

Theorem 1 (NDO Theorem): Let Zi and X 2 be two initial sets of locations 
of bubbles and let P be any program (of the type discussed above). Then 
| X p n X p 2 I ^ I Xx f\ X 2 1, where \ X \ is the cardinality of the set X. 

The NDO theorem precludes the possibility of certain types of 
computation using this model of bubble interaction. For example, it is 
impossible to have a program P such that {a, b} p = {c, d} and {a} p = 
{e}. Another consequence of the NDO theorem is the nonexistence of a 
replicating program. That is, there exists no program P* such that if 
S and S' are nonintersecting sets of vertices and X c S, P* creates a 
set X' c S' without disturbing X, so that there is a one-to-one corre- 
spondence between X and X f . In order to be able to compute all com- 
binational functions, it may be necessary to do one or both of the 
following: (i) Include other types of interactions, (ii) Code the inputs 
and the outputs. 

Let us consider the following set of instruction types, where the 
first type is the one we have discussed above. 



Name 



Notation 



Resultant Bubble Locations 



I. 



II. 



e = (a, b) 
e=(a, b) s 
e = (0, a) 
e=(l, a) 



X e = 



(X - {a})U {b} if at X,btX; 
X otherwise. 
IU \b\ if atX; 
X otherwise. 
X e = X - fa). 



X' = 



I'=IU {a}, 



Bubble 

transfer 
Bubble 

splitting 

III. Bubble 
annihilation 

IV. Bubble 

creation 

In the above set, type II instructions are necessary for replication. 
Type III instructions have the property that \X e \ < \X\ and are there- 
fore necessary to compute {a, b} p = {c}. Instructions of type IV are 
necessary for performing computations of the type (p p = {a}, where <p 
is the null set. 

Theorem 2: The four types of instructions (viz., bubble transfer, 
splitting, annihilation, and creation) are not sufficient for performing 
all computations. 



1704 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

Proof: Consider a program P such that X p = S — X. Let the first 
instruction of the program be e x . We show that a program P for comput- 
ing S — X using only the four types of instructions discussed above 
does not exist, by showing that for each type of instruction there exist 
two sets X, 5^ X 2 such that X[' = X 2 '. Since S — Xi 9± S — X a , 
the program cannot compute S — X for all X. 
Case 1: Let e x = (a, b). 

Let X 2 = (X t - {a}) *J {&}; X, , X 2 C S; a t X 1 ; b 4 X, . Then 

X\' = X 2 * = X 2 and X 1 = X 2 . 

Case 2: Let e x = (a, b) s . 

Let a e X, , b 4 X t , and X 2 = X, U {&}. Then 

Xp = X2* = X 2 and X x = X 2 . 

Case 3: Let e, = (0, a). 

Let a e Xi, and X 2 = Xi — {a}. Then 

XI' = X 2 ' = X 2 and X x = X 2 . 

Case 4'- Let e x = (1, a). 

Let a e Xi, and X 2 = X x — {a}. Then 

Xr = XS' = X 1 and X^ = X p 2 . 

Thus there exists no program P using only the four types of instruc- 
tions that can compute S — X for all X C S. This implies that comple- 
mentation cannot be done (using only these types of instructions) 
if the value of a Boolean variable is represented by the presence or 
absence of a bubble in a particular location. 

The other approach mentioned earlier for improving the computa- 
tional capabilities of magnetic bubbles is the use of suitable codes. 
Graham 4 has used two bubble locations x and x x to represent the 
value of a binary variable x. x = is represented by x — 1, Xi ■= 
and x = 1 by .To = 0, Xi = 1. However, Graham has shown that this 
representation and the bubble transfer instruction are not sufficient 
for realizing all Boolean functions. By enumerating the number of 
distinct programs, he has proven the following theorem. 
Theorem 3: There exists a Boolean function of 11 variables which 
cannot be realized by a 'program of bubble transfer instructions. 

It is not known whether there exists some coding of the variables 
which permits the realization of all Boolean functions using only the 
bubble transfer instruction. However, the four types of instructions 
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together with a coding of the type used by Graham are sufficient for 
realizing all Boolean functions. 

Theorem 4-' All combinational functions can be computed using the 
jour types of instructions and a suitable coding of inputs. 

Proof: Any combinational function can be realized using the following 
method. Let the set of all possible bubble locations be denoted by M, 
where M = SU T and S C\ T = <p. As before, values of the input vari- 
ables are represented by bubbles in subsets of S, and let each variable 
be represented by two locations .r„ and Xi . Let T be the set of possible 
bubble locations that serves as temporary storage. Any given combina- 
tional function / and its complement / are expressed in the sum of prod- 
ucts form. The following procedure is used to realize the function /: 

(i) Clear the output locations / and f, by (0, /n)(0, /i). 
(ii) Clear temporary storage by (0, 2/,o)(0, yu) for all y ia , Vn e T. 
(Hi) Copy the values of the input variables into temporary storage 
using bubble splitting instructions (x i0 , f/,- )s(#.-i , Vn)s for 
all x i0 , Xn e S. 
(iv) If x*i is a variable appearing in a product term, where x* { rep- 
resents Xi or Xi , the term J\ ?* * s computed by the set of 
instructions (y, a , y jb ) for all j 9^ i such that x*- is contained 
in the product term and a = if x* { = x, ■ , a = 1 if x* { = x, ■ , 
b — if x* t = Xi , and b = 1 if x*- = Xj . 
( v ) (Vin , /i) puts the or of all terms computed so far in /, . 
(vi) If there are additional terms, go to step ii and repeat for next 

term. 
(pit) If there are no more terms in /, repeat steps ii through vi for 
/, replacing step v by (y ia , /,,). 

The above procedure merely shows that all combinational functions 
can be computed using the four types of instructions and the par- 
ticular type of coding used. Both the program and the coding could 
be made more efficient. 

The property of the code that made the computation of all com- 
binational functions possible is the elimination of the need for com- 
plementation, which was shown to be impossible to perform with the 
four types of instructions used. A more general code which has this 
property is an m-out-oi-n (m/n) code 5 where n bubble locations, out 
of which exactly m contain bubbles, are used to represent each input 
combination. A procedure analogous to that given above may be used 
for computing any combinational function using this code. The case 
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discussed above, where each variable is represented by two bubble 
locations, is a special case of the m/n code and is referred to as the 
autosynchronous code. 5 

The m/n code is more efficient than the autosynchronous code in 
terms of the number of bubble locations required for encoding a given 
number of binary variables. Whereas 2k locations are necessary to 
represent k binary variables, the number of locations required with an 
m/n code is such that Q ^ 2*. Since m = [n/2], where [x] is the integral 
part of x, gives the largest value of O, we have ( ln " 2 i) ^ 2*. Using 
Stirling's approximation, we have [n — \ log 2 n] ^ k. For large values 
of k, the [n/2]/n code requires fewer locations than the autosynchronous 
code as the following sample values indicate: 

k ?i([n/2]/n code) 2k (autosynchronous code) 

14 16 28 

29 32 58 

61 64 122 

A given type of instruction is said to be computationally complete 
if all combinational functions can be computed using programs con- 
taining only instructions of the given type, assuming that any loca- 
tion may be initially cleared and a bubble may be initially inserted in 
any location. That is, instructions of the type (0, x) and (l,x) may be 
used for initialization, prior to the application of the program. Note 
that the bubble transfer instruction discussed earlier is not computa- 
tionally complete under our definition, although bubble transfer and 
bubble splitting together are sufficient to compute any combinational 
function. (The need for clearing locations within the program in the 
procedure discussed earlier can be eliminated by using a sufficient 
number of temporary locations.) We shall now consider some com- 
putationally complete instruction types that may be realizable by 
bubble interactions. 

Consider an instruction e = (x, y) d , defined by 

X' = X - \y\ \lx,y*X 

= {X- \x\)\J [y\ HxzX,yiX 

= X otherwise. 

If the presence or absence of a bubble in a location is treated as the 
value of a binary variable, and if the value of a variable after the 
application of an instruction (x, y) d is denoted by the corresponding 
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primed variable, the instruction can be represented by the following 
equations : 

x' = xy and y' = x © y 

where © represents exclusive-OR. 

We prove that the instruction (x, y) d is computationally complete 
by showing that we can perform duplication and also realize the nand 
function using only instructions of this type. It is well known that any 
combinational function can be realized using only two-input nand 
gates if duplication (fan-out) is allowed. Denoting locations that initially 
contain bubbles by lower-case letters with the subscript 1, the nand 
of two variables x and y can be realized by the sequence of instructions 
(x, y) d (x, aOd ; which results in, x' = xy, y' = x © y, a{ = 1 © xy = xy. 
Duplication can be performed by the sequence (x, 6i) d (6, , Ci) d , resulting 
in x' = x, b[ = 1 © x = x, c[ = 1 © x = x. Thus all combinational 
functions can be computed without special encodings. 

A timing problem associated with this instruction* must be noted. 
Two applications of the instruction (x, y) d results in xy(x © y) = 
in location / and xy © x © y = x + y in location y. Physically, if 
x = y = 1 when the instruction (x, y) d is applied, the bubble in y is 
destroyed. If the field effecting the execution of the instruction is still 
applied, the bubble in location x is transferred to y. Thus, this instruc- 
tion requires application of the field for a fixed interval of time. We 
may say that the instruction is pulse-width dependent. For each of the 
other instruction types considered, application of the same instruction 
a second time does not result in any further change in either location. 

Another type of computationally complete instruction is the condi- 
tional transfer* denoted by e = (x, y)z and defined by 

X' = {X- [x})KJ [y\ \ix,ztX,yiX 

= X otherwise. 

Duplication can be performed by the instruction (a, , a )x where a 
and a, denote an empty location and a location with a bubble respec- 
tively. After the application of the instruction, locations a and x will 
contain the variable /. Similarly, the sequence (x, y)bi , (6 S , b )x, 

* In practice, there exists an interaction involving three locations x, y, z, such 
that if y = 0, //' = x © z\ x' = z' = xz. The timing problem mentioned above 
does not occur in this case. 

t The interaction involving three locations, mentioned in the previous footnote, 
may be represented by the following sequence of conditional transfer instructions, 
if y =0 initially: (x, y)z; (z, y)x. The order of execution of these two instructions 
is immaterial. 



1708 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

yields the function xy in location b y . All combinational functions can 
be computed using only conditional bubble transfers without special 
encodings. 

III. TIME AND SPACE CONSIDERATIONS 

In the preceding section, we showed how any combinational function 
can be computed using magnetic bubble interactions, assuming that 
interactions are possible between any pair of bubble locations. Since 
interactions can occur only between bubbles in physically adjacent 
locations, this implies that bubbles that are required to interact are 
brought into adjacent locations prior to the application of the instruc- 
tion. The time taken for computing a function depends not only on the 
number of instructions in the program but also on the layout of the 
bubble locations. It is also necessary to be able to move a bubble from 
any arbitrary location to any other location without affecting the 
bubbles in other locations. We shall examine these problems in this 
section. 

3.1 Bounds on Memory Requirements 

Consider the layout shown in Fig. 1. If the locations in the shaded 
area are used as bubble locations for a program, the locations in the 
unshaded part can be used as transit points. These memory locations 
are initialized so as not to contain bubbles. Denoting the shaded and 
unshaded regions by S and T respectively, an instruction of the form 
(a, b) a, b t S will be replaced by a sequence of instructions 

(a, x 4 )(Xi , x i+l ) • • • (z ( _i , Xj)(Xj , b)(xj , Zf-t) ■ ■ ■ (x i+l , x,)(z. , a) 

where 

are adjacent locations in T so as not to affect any other location in S. 
With this layout 3n/2 — 1 locations are required for n locations which 
are useful for computation. (We shall refer to these as data points.) 
If p is the total number of locations, the maximum number of data 
points can be shown to approach 2p/3 asymptotically (K. C. Knowl- 
ton, private communication). 

The maximum path length (one-way) between two data points is 
n/2. The maximum path length can be greatly reduced at the expense 
of a slight increase in the required memory (i.e., total number of loca- 
tions) by the layout shown in Fig. 2a. The total number of locations 
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Fig. 1 — Memory layout. 

is (n/2/c - 3) -3ft + 12fc = 3n/2 + 3fc, and the number of data points 
is 2 -3ft + (ft - 1)- (re/ft - 6) + 2(n/2fc - 3) = re. The maximum 
path length (one-way) is re/2ft + 3ft - 1. [By changing the corners 
as shown in Fig. 2b, the maximum one-way path length can be reduced 
by 2, and the total number of locations can be reduced by 4 
(M. D. Mcllroy, private communication).] The maximum path length 
is dependent on ft. Since ft must be an integer and 2ft_must divide n, 
the maximum path length is minimal when ft = \/n/Q if this is an 
integer. Thus if n = 6ft 2 the maximum path length for the layout of 
Fig. 2a is (V6 \/n — 1). Some typical values of re and the minimum 
maximal path length for various values of ft are shown below: 



path length 




Fig. 2 — Memory layout to reduce maximum path length. 
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In order to obtain a lower bound on the maximum path length for 
a memory with n locations, allowing interactions between any pair 
of locations, consider the closest packing of these locations in two- 
dimensional space, which is a circle. If the radius of the circle is r, 
irr 2 = n or r = s/n/w. The maximum distance between two points 
is 2r = "vAn/V. Clearly, this bound cannot be attained because all 
points in the circle are treated as data points in deriving the bound. 
Thus the best we can ho pe t o do is to obtain a maximum path length 
of C \/n where C > \Z4/ir. The coding of Fig. 2 has maximum path 
length < s/q \/n and hence requires approximately twice as much 
time as the lower bound. 

The memory layouts discussed above allow interactions between 
all pairs of data points and are therefore suitable for realizing any 
arbitrary program. For realizing any specific function, interactions 
would be necessary only between a subset of these pairs of data points 
and consequently a smaller memory would suffice. The maximum path 
length can also be reduced in general. In this case, the memory require- 
ments and maximum path lengths derived above serve only as upper 
bounds. 

3.2 Number of Instructions for Realizing a Function 

The number of instructions and bubble locations required for realizing 
any given function using any given type(s) of instructions may vary 
over a wide range depending upon the realization, as shown in the 
following examples. 

Exam-pie 1: Consider the computation of / = x © y, using bubble 
transfer instructions and bubble splitting, if necessary. Let two loca- 
tions x and x x be used to represent each variable x, as discussed earlier. 
Using the canonical realization in the sum of products form discussed 
earlier, four minterms (two each for the function and its complement) 
have to be computed. Each minterm is computed using four bubble 
splitting instructions and one bubble transfer instruction. Finally, 
two bubble transfer instructions are required to compute the sum of 
two minterms. Thus the total number of instructions in the program 
will be 5 X 4 + 2 X 2 = 24. The total number of data points required 
is 10 (4 for x and y, 2 for /, and 4 temporary storage). Using the layout 
of Fig. 1, 14 bubble locations are required in the memory. 

The following program containing only six instructions also com- 
putes this function, if / and f\ are both initially empty: 

P = (*i , Vi) (x , y ) (x! , f ) (x Q , /„) (y l , y ) (Vi , /i). 
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Six data points are necessary, and these can be obtained using eight 
locations and the layout of Fig. 1. Since interaction between all pairs 
of data points is not required by the program, it can be shown that 
seven locations are necessary and sufficient to permit the required 
interactions. However, the following program which also contains six 
instructions and requires six data points can be realized with only six 
bubble locations: 

P' = (xi , yi)(x , y )(xi , x Q )(xo , U)(y , yi)(ya , U). 

Figure 3 shows the memory layout for this program. 

This example shows that the canonical realization may be inefficient 
in the number of instructions, computation time, and also memory 
requirements. The inefficiency of the canonical realization becomes 
even more apparent if we compare the canonical realization of / = 
Xi © x 2 © • • • © x n with the realization in the form [(x v © x 2 ) © 
x 3 ■ ■ ■ ©4 where each exclusive-OR operation is performed by a 
program of the type given above. 

The above example points out several open problems associated 
with the design of efficient programs for computing combinational 
functions. One problem is that of obtaining a program with the 
smallest number of instructions for computing a given function. For 
a given program, the assignment of locations so as to minimize the 
total number of locations required for its implementation is also an 
important problem. Since the number of memory locations required is 
dependent not only on the number of instructions, but also on the 
specific program, another open problem is that of obtaining a program 
with minimum memory requirements. 

The problem of minimizing the number of instructions in a program 
for computing a combinational function corresponds to that of finding 
a realization of the function using the smallest number of modules of 
given types. For example, the bubble transfer instruction can be repre- 
sented by the module shown in Fig. 4a. The bubble location at which 
each output appears is given in parentheses. Note that fan-out is not 



II 


I 


fo 


y, 


y 


f, 



Fig. 3 — Memory layout for program F. 
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allowed in the realization. Bubble splitting is represented by the 
module of Fig. 4b. Instructions of the type (x, y)d and the conditional 
transfer (x, y)z can be represented by modules of Fig. 4c and d 
respectively. Creation and annihilation of bubbles can be represented 
by appropriate constant inputs (0 or 1). If a realization of the func- 
tion can be obtained using a subset of these module types, then the 
realization can be readily translated into a program which computes 
the function. All bubble locations used in the program will correspond 
to inputs of the circuit. Although efficient realizations of functions 
using only modules representing the allowed bubble interactions lead 
directly to efficient bubble programs, there are at present no known 
algorithms for the former. The following example shows how a pro- 
gram is obtained from a modular realization of a function. 
Example 2: Consider the function / = ac + bed + bed, which is to 
be computed using the bubble transfer instruction only. Let each 
variable be represented by two locations. Thus, we wish to realize / 
and f using only modules of Fig. 4a, with no fan-out. The inputs to 
the circuit consist of the variables and their complements. / = 
cd + be + acd -f- dbc = c(b -\- d) + ac(b + d). We obtain the circuit 
of Fig. 5 (by trial) which realizes both / and /. Translating each module 
of Fig. 5 to the corresponding instruction, we obtain the following 
program: 

P = (b , d ) (d , c,) (6, , di) (a , c ) (d t , a„) (d Q , d t ) 

(ci , 6i) (b , a ) (c , a0 (b , c ) (d , c ). 

At the end of the program / and / will be in c and d x respectively. 
In translating the circuit to the program, no module should be trans- 
lated to the corresponding instruction until the modules connected 



IX) 



lyj 



xy 



x + y 



(x) 



(y) 



x+y 



lb) 



[cl 



(y) 



xy 



x©y 



(X) 

(z) 



xy + xz 

xz+y 

z 



id; 



Fig. 4 — Modules representing instructions: (a) Bubble transfer, (x, y). (b) 
Bubble splitting, (x, y)a. (c) Type (x, y)t. (d) Conditional transfer (x, y)z. 
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Fig. 5 — Circuit for example 2. 

to its input terminals have been translated. It is also important to 
observe the proper correspondence between bubble locations and the 
output terminals of a module. After applying an instruction (x, y), 
xy appears at location x, and x + y appears at location y. 

The other types of modules shown in Fig. 4 can be used in a sim- 
ilar manner for deriving programs using other types of instructions. 
However, there are no known algorithms for realizing any arbitrary 
function using only modules of a given type, so as to minimize the 
number of modules used. When such algorithms become available they 
can also be used for obtaining efficient programs of magnetic bubble 
interactions. 



3.3 Speed-up of Computations 

We have seen how a combinational function can be computed by 
different programs having different numbers of instructions and 
memory requirements. If the instructions in a program are executed 
one at a time, as was implicitly assumed so far, then the time for 
computing a function depends on the number of instructions in the 
program. However, it may be possible to speed up the computation 
by executing several instructions simultaneously. 

An obvious method of obtaining parellelism in a program is by 
executing independent instructions simultaneously. Instructions that 
may be executed simultaneously can be determined from the circuit 
corresponding to the program, discussed in the preceding section. The 
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modules in the circuit are arranged according to levels as follows: 
Modules whose inputs are only input variables (or their complements) 
are assigned to level 1. The level of any module is defined to be the 
smallest integer greater than the levels of the modules connected to its 
input terminals. With levels assigned to all modules in the circuit in 
this manner, the instructions corresponding to these modules are 
executed in the order of the level numbers, and the instructions cor- 
responding to all modules in one level may be executed simultaneously. 
However this may necessitate a further restriction on a valid memory 
allocation scheme. 

Example 3: Consider the realization of the function / = ac + bed + bed 
in Fig. 5. The modules can be arranged according to levels as follows: 
level 1-(1, 3, 4); level 2-(2, 5, 9); level 3-(6, 7, 8); level 4-(10); level 
5— (11). The program P can be executed in five steps as follows: 

(i) (b , d ) (bi , dj) (a , c ). 

(ii) (d , c s ) (di , a ) (c , Oi). 

(Hi) (do , di) (d , bx) (b , a ). 

(iv) (b , Co). 

(v) (c, , c ). 

Another situation where simultaneous execution is possible becomes 
obvious by considering the instructions (a, b)(a, c). If these instructions 
are executed simultaneously, the location a will contain dbc, whereas 
the contents of locations b and c are indeterminate because they depend 
on which of the two instructions is actually executed first. However, 
if we are interested only in the product term abc, the two instructions 
may be executed simultaneously. In general, all (n — 1) instructions 
for forming a product of n variables can be executed simultaneously. 
Similarly, the simultaneous execution of the instructions (a, s) (b, s) ■ ■ ■ , 
leaves the sum a + b + • • ■ + s in location s with the contents of 
a, b, • • • indeterminate. Using this technique, the program of Example 3 
can be executed in four steps by replacing the instructions of steps 
iv and v by (b , c Q )(ci , c ), executed simultaneously. Instructions of 
the type (x, y) d can also be executed simultaneously with similar results. 
For instance, the simultaneous execution of (a, b) d and (a, c) d leaves 
the product abc in location a, with the contents of locations b and c 
indeterminate. Similarly, if (a, c) d and (b, c) d are executed simul- 
taneously, the location c will contain a b c, but the contents of 
locations a and b will be indeterminate. 
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It may also be possible to speed up computations in a manner sim- 
ilar to the use of completion signals in asynchronous circuits. 5 " 7 A 
combinational function / can be expressed as a sum of disjoint 
minterms. The function can be computed by computing these min- 
terms, one at a time. As soon as a 1 is generated, the computation is, 
in effect, complete because the value of the function is known to be 1. 
The computation can be terminated if the presence of a bubble in 
the output location can be detected by the control circuitry. The con- 
ditional transfer instruction (.r, y)z discussed earlier may be useful 
for this purpose. 

3.4 Iterative Array Realizations 

Our discussions so far have been restricted to realizations which 
minimize the number of instructions or the computation time. In this 
section, we shall consider some techniques for obtaining regular struc- 
tures, which are capable of performing computations by the applica- 
tion of uniform magnetic fields. In such a structure, the same type 
of instruction is executed simultaneously in all parts by the applica- 
tion of a uniform magnetic field, the type of instruction being con- 
trolled possibly by the direction of the magnetic field. 

Since regularity of interconnections is a highly desirable feature in 
integrated circuits also, a great deal of effort has been directed toward 
the realization of functions as iterative arrays. 8 - Such arrays may 
be classified as uniform cell function arrays (in which the function 
realized by each cell is identical and different functions are realized 
by the array by changing the inputs to some cells) and nonuniform 
cell function arrays in which different cells realize different functions, 
depending on the position of the cell in the array. The operation of 
uniform cell function arrays can be simulated by magnetic bubble 
interactions with uniform magnetic fields. 

The rectangular array shown in Fig. 6a may be used for realizing 
any combinational function of n variables by merely changing the 
connections to the last row of cells. A typical cell in the array is 
shown in Fig. 6b. It can be shown that the vertical outputs of the n th 
row consist of the 2" minterms of n variables. Since at most one of 
the minterms will be 1 at any time, the last row will compute the 
sum (or) of all minterms connected to it. Any function can be realized 
by connecting to the cells in the last row only those vertical outputs 
of the n th row that correspond to 1-points of the function. 

Note that the number of cells in the array is (n + 1)2". The longest 
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Fig. 6 — (a) Rectangular array, (b) Typical cell in array. 

path in the array is 2" + n + 1. Therefore, the time required for 
computing a function of n variables is proportional to 2" + n + 1. 

The operation of this array can be simulated by magnetic bubble 
interactions as follows: Let each cell in the array be replaced by two 
bubble locations as shown in Fig. 7. Let the application of a magnetic 
field to a cell cause the execution of the instruction (x, y) d and transfer 
of the resultant bubbles as shown by the arrows in Fig. 7. (It may be 
necessary to apply a magnetic field in one direction followed by a field 
in a different direction to accomplish this.) Let the cells in the ith row 
and jth column be denoted by C,, . If the left and right locations in 
Ca represent x and y respectively, application of the field to C,,- results 
in xy in the right location of C I+1 ,,- and x © y in the left location of 
Ci,j+i . Initially, we set the left locations of all cells in the first column, 
except the cell in the last row, to correspond to the values of the n 
variables. The left location of the first cell in the last row should be 
empty. The right locations of all cells in the first row initially contain 
bubbles. Paths from the nth row to the (n + l)st row which correspond 
to minterms for which the function is are inhibited. At t = 1, the 
field is applied to Cu . At t = 2 the field is applied to Ci 2 and C 2 i ; 
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at t = 3 to C, 3 , C 2 2 , and C 3l ; and so on. That is, at t = k, 1 ^ A; ^ 
2" -f- ra, the field is applied to all cells C, such that * + j = fc + 1. At 
£ = 2" + fl. + 1, the value of the function will be in the right location 

of C n+1 , 2 » • 

Instead of applying fields to successive diagonals, the fields may be 
applied to the entire platelet. Now, the right locations of all cells in 
the top row should be connected to "bubble generators," so that they 
always contain bubbles. The values of the n variables, contained in 
the left locations in the first column, may be changed every n units 
of time. The output corresponding to any input combination will 
appear at the output location 2" + n + 1 units of time after the 
application of the input combination. If input changes occur n units 
of time apart, the correct outputs will also appear n units of time 
apart. However, the contents of the output location will not be correct 
between these fixed instants of time. Therefore, it is necessary to read 
the output at the appropriate instants of time. 

The preceding discussion shows how a uniform cell function iterative 
array can be simulated by magnetic bubbles by applying identical 
sequences of instructions to successive sets of bubble locations. This 
provides us with a systematic way of realizing all combinational func- 
tions. If conditional operations are permitted, nonuniform cell func- 
tion arrays can also be simulated. Such arrays can be made smaller 
than uniform cell function arrays. The number of cells required for 
realizing any arbitrary function of ?i variables is proportional to 2 n 
with nonuniform arrays compared to n-2 n required for uniform cell 
arrays. 8 However, the cells of the nonuniform array are likely to be 
more complex than those in a uniform array. 

It is interesting to compare the size of array realizations to the size 
of realizations without any restrictions on the interconnection struc- 
ture. Bounds on the latter will provide us with an indication of the 
space and time requirements for computing an arbitrary combina- 
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Fig. 7 — Simulation of array by magnetic bubble interactions. 
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tional function with magnetic bubbles. D. E. Muller 10 has shown that 
for any set of elements capable of realizing all combinational func- 
tions, the number of elements N for realizing any arbitrary function 
of n variables satisfies the inequality 

^-" < N < ^ , 
n n 

where C\ and C 2 are constants whose values depend on the set of 
elements used. The proof of the upper bound is constructive and uti- 
lizes only one type of cell. It is outlined below because of its possible 
applicability to magnetic bubble computations. 

The cell used realizes a function of three inputs ab -f- dc. Any function 
of k variables can be expressed as 

}\Xi , X 2 , • • • , Xk) = Xk]\Xi , X2 , " " ■ , Xk-\ , 1) 

+ x k j(xi , x 2 , •■• , x k -i , 0). 

If all functions of k — 1 variables are available, every function of k 
variables can be realized using exactly one cell. There are 2 2 functions 
of k variables (which also include all functions of k — 1 variables, 
k — 2 variables, etc). Thus 2 2 cells are sufficient for realizing all func- 
tions of k variables. 

By expanding about any variable, the output function can be formed 
from two functions of n — 1 variables. Repeating the procedure, it 
may be formed from 2 n_ * functions of k variables, using 2 n ~* — 1 ele- 
ments. Since 2 2 ' elements are sufficient for realizing all functions of 
k variables, 

N = 2 2 ' + T~ k - 1. 

The value of k may now be chosen so as to minimize N. An approximate 
minimum is N < C 2 (2"/n), for some constant C 2 . Thus the bound 
on the number of cells required is greatly reduced for nonregular 
structures. 

IV. SUMMARY 

Three different mathematical models of magnetic bubble interac- 
tions were studied and each of them was shown to be sufficient for 
computing all combinational functions. Some methods of speeding up 
computations were presented. Geometrical patterns in which it is 
possible to move bubbles between any pair of data locations and also 
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minimize the maximum path length were examined. A uniform struc- 
ture for computing arbitrary functions was also presented. These 
structures have the advantage that computations can be carried out 
by the application of uniform magnetic fields to the entire platelet. 
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